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Abstract. Results related to the possible chemical equilibration of hadrons in heavy 
ion collisions are reviewed. Overall the evidence is very strong with a few clear and well- 
documented deviations, especially concerning multi-strange hadrons. Two effects are 
considered in some detail. Firstly, the neglect of (possibly an infinite number) of heavy 
resonances is investigated with the help of the Hagedorn model. Secondly, possible 
deviations from the standard statistical distributions are investigated by considering 
in detail results obtained using the Tsallis distribution. 
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After analysing particle multiplicities for two decades a remarkably simple picture has 
emerged for the chemical freeze-out parameters [Il[2l[3]. Despite much initial skepticism, 
the thermal model has emerged as a reliable guide for particle multiplicities in heavy ion 
collisions at all collision energies. Some of the results, including analyses from [U El El [7] , 
are summarised in Fig. [1] Most of the points in Fig. [1] (except obviously the ones at 
RHIC) refer to integrated (Att) yields. A clear discrepancy exists in the lower AGS 
beam energy region between the (published) mid-rapidity yields and estimates of the 
Att yields. The latter tend to give higher values for the chemical freeze-out temperature. 
This will have to be resolved by future experiments at e.g. NICA and FAIR. When 
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Figure 1. Values of the freeze-out parameters obtained at beam energies ranging from 
1 GeV to 200 GeV 
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the temperature and baryon chemical potential are translated to net baryon and energy 
densities, a different, but equivalent, picture emerges shown in Fig. O This clearly 
shows the importance in going to the beam energy region of around 8-12 GeV as 
this corresponds to the highest freeze-out baryonic density and to a rapid change in 
thermodynamic parameters [H [9] . 

The dependence of /i^ on the invariant beam energy, ^/snn, can be parameterized 
as [3] 

_ 1.308 GeV 

^^^^^ - 1 + 0.273 GeV-^^- 

Similar dependences have been obtained by other groups |Tl|2]. and are consistent with 
the above. This predicts at LHC /i^ ~ 1 MeV. 

To analyze the changes around 10 GeV use can be made of the entropy density, 
s, divided by which has been shown to reproduce the freeze-out curve [3J very well. 
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Figure 2. The hadronic freeze-out hne in the Pb — ^* phase plane as obtained from 
the values of /is and T that have been extracted from the experimental data in [3]. 
The calculation employs values of /ig and /is that ensure (S) = and (Q) = 0.4(B) 
for each value of /i^. Also indicated are the beam energies (in GeV/N) for which the 
particular freeze-out conditions are expected at either RHIC or FAIR or NICA. 
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Figure 3. Values of entropy density divided by following the chemical freeze-out 
values. 



This allows for a separation into baryonic and mesonic components, shown in Fig. [3l it 
can be seen that mesons dominate the chemical freeze-out from about ^/sNN ~ 10 GeV 
onwards. 
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Figure 4. Regions in tlie T—^b plane where baryons or mesons dominate as indicated. 



2. Resonance Gas and Hagedorn Spectrum 

It is possible to compare analytically the resonance gas which uses a finite number of 
resonances up to some maximum mass, with a Hagedorn gas which contains an infinite 
number of resonances with an exponentially rising number of resonances as the mass 
increases. It is well known that at some point the Hagedorn resonance gas will show a 
divergence when the temperature reaches the Hagedorn value [HI [121 [13 [H] ■ 
The speed of sound is given by 

s/n 

where the derivative is taken at a fixed entropy per particle. 

For an ideal Boltzmann gas of identical scalar particles of mass tuq and three charge 
states ( "pions" ) contained in a volume V, the grand partition function is defined as 

-1 N 

- J d^p exp{-y^p2 _^ rno^/T} 



3 V 



_(27r)3. 

This expression can be evaluated, giving 

VTm ^ 

lnZ(T,y) = 3^-— ^ir2(mo/T), 

The speed of sound in this gas is given by 

i_3 ^ mo^ir2(mo/T) 

Now extend this to an ideal Boltzmann gas of resonances [IT], described by an 
exponentially increasing mass spectrum of the Hagedorn form 

p(m) = 35(m — mo) + Am~^ exp{m/Tc} 9{m — 2mQ). (2) 




Figure 5. Speed of sound calculated using the thermal model, following the values 
of the chemical freeze-out curve but with different contributions to the resonance gas 
determined by the masses of the resonances. 



The ^-function assures that the resonance spectrum starts above the two-pion threshold. 
The grand partition function is now given by 

VT 

In Z{T,V) = —SmlK^imo/T) (3) 
VT r°° 

H -A / dm exp{m/Tc} K2{m/T), (4) 

271^ J 2mo 

and the speed of sound by 

1 3mf,K2(mo/T) + A [°° dm exp{m/Tc} Kiim/T) , ^ 

^2 + A Tj^^^dm m.-^ exp{m/Tc} [4TK2{m/T) + mKi{m/T)]' ^ ' 

The second term in the numerator diverges as T — )■ Tc, which in turn causes the speed 
of sound to vanish there [HI [121 IE]- At low temperatures there is almost no difference 
between a Hagedorn gas and a thermal model containing only a limited number of 
resonances, however at higher temperatures the calculated speeds of sound become very 
different as shown in Fig. [2J It is thus always necessary to check if results obtained in 
the thermal model are stable against the addition of a Hagedorn-type of mass spectrum. 
Fortunately, for many quantities of interest the answer is yes [T^ . 

3. Non-extensive Tsallis statistics 

For the Tsallis distribution [15], one replaces the standard expression for the entropy 
based on 

S = -Y,PMPii (6) 
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Figure 6. Comparison between tlie Boltzmann and Tsallis distributions. 



with 



^7 



(7) 



This introduces a new variable q, often referred to as the Tsalhs parameter. In the hmit 
where this parameter goes to 1 one recovers the Boltzmann entropy 



lim St 



S 



(8) 



The physical interpretation of the Tsallis parameter is not obvious, we will subscribe 
here to the one presented in Ref. [23] . 

We have repeated the analysis for particle yields using the Tsallis distribution. The 
particle densities of particle yields were calculated using [T8] : 

i^_^n-V('?-i) 



n'{E) 



T 



(9) 



A very interesting application of this distribution to the transverse momentum 
distribution observed in heavy ion collisions has been presented at this conference by 
the STAR collaboration [20]. In the limit where the parameter q tends to 1 one recovers 
the Boltzmann distribution: 

E - fi' 



lim (E) 



exp 



T 



(10) 



A comparison between the two distributions is shown in Fig. [3] Clearly, at some value 
of q, an integral over a Tsallis distribution will no longer give a convergent result. 

A possible interpretation of the Tsallis Parameter q has been presented in [23]. One 
starts by rewriting the Tsallis distribution as a superposition of Boltzmann distributions 
with different temperatures, this is possible using a distribution function g: 

'E - "\\ -1/(9-1) 



f{E)= l + (g-l) 



T 



:ii) 
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Figure 7. The x^/d.o.f. of the fits as a function of the Tsallis parameter q. 



The precise form of the function g has been given in Ref . [23] . The average temperature 
is then given by the T parameter appearing in the Tsalhs distribution: 




= ^ (13) 
and the Tsalhs parameter q is the deviation around this average Boltzmann temperature, 




Thus in the hmit where q goes to one, this goes to zero. 

The resuhing value of shows an interesting dependence on the parameter g, as 
shown in Fig. |3l It must be added that most of the thermodynamic parameters show a 
very strong dependence on the parameter q which necessitates a complete review of the 
physical picture behind chemical freeze-out. The temperature is shown in Fig. [HI Other 
variables like the volume and the chemical potential also show a strong variation with 
the Tsallis parameter q [T8] . 
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Figure 8. The freeze-out temperature as a function of the Tsahis parameter q. 
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